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We present some results for SU(2) with one adjoint Dirac flavour from lattice studies. 
Data for the spectroscopy, the static potential, topological charge, and the anomalous 
dimension of the fermionic condensate are included. Our findings are found to be in¬ 
consistent with conventional confining behaviour, instead pointing tentatively towards a 
theory lying within or very near the onset of the conformal window, with an anomalous 
dimension of the fermionic condensate of almost 1. Implications of these findings on 
the building of models of strongly-interacting dynamics beyond the standard model are 
discussed. 


1. Introduction 

Novel strong interactions provide a promising explanation for the mechanism of 
electroweak symmetry breaking. Theories near the onset of the conformal window 
with an anomalous dimension of the chiral condensate of order one may in particular 
be relevant for phenomenology. It is important that these properties are verified for 
candidate models in a first-principles, fully non-perturbative fashion, and numerical 
lattice simulations provide an indispensable tool to this end. Recent lattice calcu¬ 
lation's^ have demonstrated that the SU(2) gauge theory with two adjoint Dirac 
flavours is in the conformal window, while studies of many-flavour qCDIMI and of 
the SU(3) gauge theory with two Dirac flavours in the sextet representation^StH are 
ongoing. 

In these proceedings we present results from an ongoing investigation of the 
SU(2) Yang-Mills theory with one adjoint Dirac flavour, using Monte Carlo studies 
of the theory discretised on a spacetime lattice. The chiral symmetry breaking of 
this model produces only two Goldstone bosons, making it insufficient to give mass 
to the and Z bosons. Thus the theory needs to be extended for phenomeno¬ 
logical considerations. The main purposes of this study are firstly, to offer potential 
insights into the connection between near-conformal gauge theories and the pres¬ 
ence of a large anomalous dimension. Secondly, to seek to pin down the lower end 
of the conformal window, which may give insight as to whether an interpolating 
theory between one and two adjoint Dirac flavours could be of phenomenological 
interest. Third, as an end rather than a means, as investigation and classification 
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of gauge theories is also of theoretical interest. 


2. Lattice computations 

We simulate the theory on a 4-climensional hypercubic lattice, using the Wilson 
gauge action and the Wilson fermionic action, giving 

P x,y 

where the first sum is over plaquettes, and D(x,y) is the massive Dirac operator: 

y) = Sx,y — K [(1 — Ufj,{x)6y^x+fi + (1 + Jy.) Uj^{x — fi)dy^x—fi\ ■ 

K parameterises the bare fermion mass am, being related to it as k = 1/(8 + 2am). 
(3 = 2N/g‘^ = 4/g^ is the inverse coupling, where TV = 2 is the dimension of the 
gauge group. 

Gauge ensembles were generated using the HiRep code suite, making use of 
the Rational Hybrid Monte Carlo (RHMC) algorithm. Further detail of the initial 
set of ensembles are given in^ the newer ensembles introduced in this work were 
generated using the same procedures. 

Ensembles were generated at two values of /3: at /3 = 2.05, masses of —1.524 < 
am < —1.475 were used, while at /3 = 2.2 the range —1.378 < am < —1.280. With 
the lattice volume fixed to Nt x iV£, spatial lattice sizes of 12, 16, and 24 were 
chosen, in all cases being cautious to avoid the effects of finite volume artefacts in 
our observables. 


Table 1. Parameter sets considered at ^ = 2.05 and 2.2, showing their volumes, bare masses, 
PCAC masses, and number of configurations. 




(a) /3 

= 2.05 




(b)/3 

= 2.2 



Nl 

—am 

(UTlPCAC 

A^conf 


Nl 

—am 

O’PT'PCAC 

A^conf 

B1 

12 

1.475 

0.1493(5) 

1500 

2B1 

12 

1.280 

0.2017(6) 

4000 

B2 

12 

1.500 

0.1113(8) 

1500 

2B2 

12 

1.290 

0.1891(5) 

4000 

B3 

12 

1.510 

0.09226(92) 

4000 

2B3 

12 

1.300 

0.1761(7) 

4000 

Cl 

16 

1.475 

0.1485(4) 

1500 

2B4 

12 

1.310 

0.1633(7) 

4000 

C2 

16 

1.490 

0.1279(2) 

1500 

2C1 

16 

1.320 

0.1487(3) 

4000 

C3 

16 

1.510 

0.09111(31) 

1500 

2C2 

16 

1.340 

0.1189(4) 

4000 

C4 

16 

1.510 

0.09048(52) 

1500 

2C3 

16 

1.350 

0.1021(5) 

4000 

C5 

16 

1.514 

0.08223(34) 

1500 

2C4 

16 

1.360 

0.0844(5) 

4000 

C6 

16 

1.519 

0.06587(37) 

1500 

2D1 

24 

1.360 

0.0830(10) 

1721 

D1 

24 

1.510 

0.09130(27) 

1534 

2D2 

24 

1.370 

0.0652(2) 

4000 

D2 

24 

1.523 

0.04722(43) 

2168 

2D3 

24 

1.378 

0.0486(5) 

1505 

D3 

24 

1.524 

0.04081(71) 

800* 







Note: * Run D3 used configurations saved every 5 trajectories rather than every one, giving 
equivalent statistics to 4000 trajectories. 


Our aim in studying this theory is to probe its infrared regime—in particular, 
we would like to identify conformal properties and calculate the chiral condensate 
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anomalous dimension. To this end, followin^^^, we have calculated the mass spec¬ 
trum in the chiral regime, and tested finite-size scaling predictions. Additionally, 
we have performed fits of the Dirac mode number^. 

The massless fermionic action admits a global SU(2) chiral symmetry, which 
is broken by a non-zero chiral condensate to SO(2) = U(l). The generator of 
the unbroken subgroup is identified with the baryon number B. Since parity is 
unbroken, the states of the theory are labeled by , where P is the parity, and by 
their spin. The spectrum comprises mesons, baryons, glueballs, and glue-fermion 
composite states^. In addition to these states, we also measure the string tension. 

Mesonic states are accessed by the usual correlation functions of operators com¬ 
posed of fermion fields, which necessarily contain disconnected contributions due 
to the single flavour. Baryonic states meanwhile are accessed via the connected 
correlation functions of appropriately-chosen operators of the same form as for the 
mesons. Mesons have baryon number B = 0, while baryons have B = ±2. 

Meanwhile glueball operators consist of linear combinations of ordered prod¬ 
ucts of link matrices on closed path, chosen to transform irreducibly under spin 
transformations, parity, and charge conjugation. 

The string tension is accessed via two methods: via operators of blocked smeared 
Polyakov loops, fitted to an effective string theory prediction, and via the tempo¬ 
ral asymptotic behaviour of expectation values of Wilson loops, giving the static 
potential, which can be fitted to a Cornell-like form. 

For further details of the methods used to generate the ensembles and to extract 
these observables, we defer to the more detailed presentation given in^^Sl. 



Fig. 1. Sample topological charge histories for ensembles (a) Cl, and (b) D2. 


As a check of ergodicity, we look at the topological charge history, illustrative 
examples of which are shown in Fig. Good ergodic movement is seen, with no 
sign of frozen regions. To verify that we are not in a region of lattice artefacts, the 
histogram of the average Polyakov loop in each direction is checked; in Fig. [^we 
show illustrative examples showing an unbroken centre symmetry via the unimodal 
distribution—a broken centre would be indicated by a bimodal distribution. 
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Fig. 2. Polyakov loop distribution for ensembles (a) C4, and (b) D2. 


3. Spectrum 
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(a) (b) 

Fig. 3. Bare spectral quantities for (a) 0 = 2.05, and (b) 0 = 2.2, as a function of ampcAC- 




0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0 0.05 0.1 0.15 0.2 0.25 0.3 

“ ”»PCAC “ ”»PCAC 

(a) (b) 


Fig. 4. Spectral quantities for (a) 0 = 2.05 and (b) 0 = 2.2, scaled by the string tension as a 
function of the PCAC mass ampcAC- all cases the behaviour observed is consistent with flat. 


We now turn our attention to the spectrum of the theory. As seen in Fig. all 
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spectral observables of the theory, including meson, baryon, and glueball masses, 
and the string tension in lattice units, decrease monotonically towards zero as 
OTOPCAC -t 0. 

In Fig. we show the spectrum of the theory, now as dimensionless ratios to the 
string tension, both for f3 = 2.05 and with available preliminary data for /? = 2.2. 
All states show constant mass ratios to two standard deviations throughout the 
range in which they were observed. At /3 = 2.05, the 0+ scalar meson is seen to be 
the lightest state in the scaling region, and appears to be degenerate with the O'*'"*' 
glueball, with which it is anticipated to mix. 


4. Chiral condensate anomalous dimension 



(7* = 0.9) ( 7 * = 1.0) ( 7 . = 1.1) 

Fig. 5. Plots of LrriQ- as a function of for /5 = 2.05 for three values of 7 *. The 

results appear to identify a universal curve for 0.9 < 7 * < 1.0. 


We consider two independent methods to extract the anomalous dimension 7 *. 
The first gives an approximate value via finite size scaling predictions. For a con¬ 
formal theory, a spectral quantity mx of the system in a finite volume of spatial 
extent L, as L —>■ 00 and the combination is kept constant, obeys the 

relation 

Lmx = / ’ 

for some unknown function /. Thus by plotting Lmx as a function of 

for our data at various volumes, we may extract 7 * by tuning it until the data lie 

on a single universal curve. 

A sample of such plots are shown in figure from this we can see a best fit in 
the range 0.9 < 7 » < 1.0. 

The second route we use to obtaining 7 * is by fitting the Dirac mode number 
Jy(fi) as a function of the Dirac eigenvalues 17^^. Since at low fl the distribution 
is distorted by the finite fermion mass, while at large D the scaling region ends, 
we look for an ensemble with a large intermediary region to examine; thus a large 
lattice and small mass makes sense. We therefore examine the D2 ensemble, and 
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Upper end of window 


Upper end of window 


Lower end of window: 


0.100309 (-+-( 0.113707 0.128894 h 

0.104589 0.118559 h-B-h 0.134394 ^ 

0.109053 0.12:1618 -:--h 0.140130 h 


H 0.146110 i-v-i 0.165625 
^ 0.152345 I Y I 0.172693 
H 0.158846 i-A—I 0.180063 




H 0.187747 i-^-( 0.212824 h 
. 0.195759 0.221906 h 


H 0.241250 h-A-H 0.273473 h-O-h 0.310000 
H 0.251546 h-V-H 0.285144 h-4-h 0.323229 


Fig. 6 . Plots of the fitted parameters A, uq, 7 * and as a function of the window bounds. A 
region of stability is identified around 7 * = 0.92(1). 


look to fit 

a~‘^T'{n) « a~‘^T'o(m) + A [(afl)^ — (aM)^] . 

Since M here is a renormalised mass and not known a •priori, we have four param¬ 
eters to fit: Vq, a, (aM)^, and 7 *, However, since the scaling region is also not 
known a priori, an additional fitting degree of freedom is introduced in the choice 
of fitting window. Since the fit is highly sensitive to initial conditions, each possi¬ 
ble window is fitted many times for varying starting parameters, and uncertainties 
estimated using the bootstrap method. We then look for a region where the fit is 
stable for a range of window lengths and positions, as can be seen in Fig. we see 
this gives a best ht of the anomalous dimension of 0.92(1). Further details of the 
fitting procedure are given in^. 

5. Conclusions 

We have extended our previous investigation of SU(2) with one adjoint Dirac flavour 
to a second value of j3, where we have a preliminary set of observables available. 
At both values of /3 we see behaviour consistent with conformality and inconsistent 
with QCD-like confinement: namely, monotonic change of all spectral quantities, 
whose ratios are flat, consistent with hyperscaling behaviour. At one value of /3, we 
find a large anomalous dimension of the chiral condensate, 0.9 < 7 * < 0.95, with a 
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best fit of 0.92(1). Work is ongoing to extend the study at /? = 2.2 to the full set of 
observables studied for (5 = 2.05—including the anomalous dimension of the chiral 
condensate, which we currently have insufficient data to determine conclusively— 
after which we intend to study sufficiently many values of j3 to allow a continuum 
limit extrapolation. 
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